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Abstract. We have analysed the optical luminosity-redshift 
distribution properties of bright QSOs, using a new large sam- 
ple from the Hamburg/ESO survey. The sample provides in- 
sight into the hitherto poorly sampled bright tail of the lumi- 
nosity function, allowing to monitor its evolution with redshift 
up to z « 3. The slope increases significantly towards higher 
z, inconsistent with the predictions of pure luminosity evolu- 
tion, but also with other recently proposed parameterisations. 
This phenomenon is opposite to what would be expected from 
gravitational lensing, showing that magnification bias does not 
significantly distort the QSO luminosity function within the 
redshift range covered. The space density of high-luminosity 
QSOs continues to increase up to the high-redshift limit of the 
survey, without indication of reduced evolution above z — 2. 
The sample also permits an improved estimate of the local 
(z ~ 0) luminosity function of QSOs and bright Seyfert 1 nu- 
clei, over the luminosity range — 27 ^ Mbj ^ —20. No evidence 
for a break or change of slope is found down to absolute mag- 
nitudes Mbj ~ -20. 



Key words: Quasars: general - Galaxies: Seyfert - Galaxies: 
evolution - Galaxies: luminosity function - Cosmology: obser- 
vations - Gravitational lensing 

1. Introduction 

The wide range of quasar luminosities, spanning altogether 
over four decades from low-activity Seyfert galaxies up to the 
most luminous objects in the universe found at high z, al- 
lows individual surveys to access only small portions of the 
luminosity-redshift-plane at given redshift. While at interme- 
diate flux levels there are now many samples available, large 
uncertainties still exist at the faint and bright tails of the QSO 
luminosity function (QLF). The lack of good constraints for the 
bright end is particularly grave in the optical waveband, where 
the Pa lomar -Green Bright Quasar Sample (BQS; Schmidt & 
Green 1983 1 has been the only significant contributor for now 
15 years, although ist is now well k nown to be substantially 
incomplete (Goldschmidt et al. |1992| ; Kohler et al. |1997| ). 

The most widely quoted analysis of the quasar luminos- 
ity function is due to Boyle and collaborators (e.g., Boyle et 



al. 1988). who concluded that the evolution proceeds along in- 



creasing luminosity, with the simple scheme of 'pure luminos- 
ity evolution' (PLE) providing a statistically acceptable fit to 
their data. However, the Boyle et al. analysis had to rely on 
the BQS and was therefore affected by the incompleteness of 
that sample. Newer surveys sensitive to the brighter part of 
the QSO population have questioned the validity of the PLE 
pictur e. In a preliminary analysis of the LBQS, Hewett et al. 
(1993) showed that the PLE model predicts too few bright and 
too many faint low-redshift quasar s. A si milar conclusion was 
reached by Goldschmidt & Miller ( 1998[ ) using the Edinburgh 
survey, by La Franca & Cristiani (1997) from a com bined anal- 
ysis of EQS and HBQS, and by Kohler et al. (1997) comparing 
their 'local' QLF with the PLE prediction. 

For redshifts z > 2, the general view is that the evolution 
slows down considerably around z — 2 and eventually reverses. 
The actual onset of this reversal is, however, not well con- 
strained, and ma y well depend on luminosity. In the PLE model 
of Boyle et al. (1991), comoving space densities are roughly 
constant for redshifts z <; 1.9, while the results of Hewett et 



al. (1993 1 and Warren et al. ( 1994) suggest that space densities 
still increase up to z — 3, although at a reduced rate, but drop 
precipitously thereafter. A marked decline in space densities to- 



wards higher z was also found by Schmidt et al. ( 1995 ). How- 



ever, all these surveys sampled mainly intermediate-luminosity 
QSOs with Mb — —26. For highest luminosities, Mb ^ —28, 
there is very little material available, and the few results so far 



are controversial. Irwin et al. ( 1991 ) claimed that there is no ev- 
idence for evolution at all, up to z = 4.5. This has been disputed 



by Kennefick et al. (1995), who obtain an order of magnitude 
lower space densities at z — 4.3 than Irwin et al. Resolving this 
controversy will require new wide-angle survey material; the 
first r esults from the just started Sloan Digital Sky Survey (Fan 
et al. 1999) indicate that data on high-luminosity QSOs in the 



very early universe will soon be available in substantial num- 
bers. 

In this paper we present the contribution of a new large 
sample of bright QSOs, drawn from the Hamburg/ESO survey 
(HES; Wisotzki et al. 119961, Reimers et al. 119961). This sam- 



ple contains many of the most luminous quasars known and 
is suited to study the evolutionary properties at the bright part 
of the QLF over a wide range of redshifts. It is also the first 
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Fig.l. Distribution of the 415 QSOs of the fiux-Hmited HES 
sample over absolute magnitudes and redshifts. 



quasar sample large enough to effectively replace the BQS, 
with higher completeness, photometric accuracy, and less af- 
fected by redshift-dependent selection biases. The main aim of 
this paper is to use the HES sample to derive constraints on 
the evolution of luminous QSOs. We do not attempt to com- 
pute new evolutionary models, as that would require merging 
the HES with other, fainter samples - a complication that we 
wish to avoid in the present paper. 

2. Observational data 

The Hamburg/ESO survey was initiated in 1990 as an ESO key 
programme, to perform a wide-angle search for bright QSOs in 
the southern sky. The survey uses objective prism plates taken 
with the ESO Schmidt telescope and digitised with the Ham- 
burg PDS microdensitometer (for more details on the survey 
see Reimers & Wisotzki 1997; Wisotzki et al. 1999 - hereafter 
Paper 1). 



The selection of QSO candidates from the database of dig- 
ital spectra involves a multitude of selection criteria. Extensive 
follow-up spectroscopy conducted at ESO has allowed to con- 
struct a new flux-limited sample of 415 objects within the red- 
shift range ^ z < 3.2, with optical magnitudes Bj ^ 17.5; 
a full discussion of the sample properties is given in Paper 1, 
of which only some essentials are summarised here. The sam- 
ple was compiled from processing 207 HES Schmidt fields 
with complete follow-up spectroscopic identification down to 
well-defined flux limits. Each field has its own limiting magni- 
tude depending on the quality of photographic plates and see- 
ing. While one single Schmidt plate formally subtends over 
~ 5° X 5°, the effective survey area is reduced because of 
overlapping adjacent plates and loss of processable area. The 
total effective area ilgff is 3700 deg^ for Bj < 14.5, and a 
monotonously falling function of magnitude for fainter objects. 
The uncertainties of optical photometry are generally smaller 
than 0.2 mag. All redshifts were determined from follow-up slit 



spectra, with typical continuum S/N ratios of ^ 20, and for no 
object the assignment of an appropriate redshift was in doubt. 

Due to its wide range of selection criteria and highly au- 
tomated surveying procedure, the HES is less affected by 
redshift-dependent selection biases than many other optical 
surveys. In Paper 1 we give an extensive discussion of the com- 
pleteness for this sample; the central points are summarised as 
follows: 

(1) The redshift distribution does not show evidence for sig- 
nificantly enhanced or reduced selection efficiency in certain 
redshift regions. 

(2) More than 99 % of the previously known QSOs located 
within the survey area and above the flux limits were recov- 
ered. 

(3) QSOs within the usual range of spectral properties are al- 
ways selected. 

(4) The surface density of bright QSOs as measured in the HES 
is higher by a factor of 1 .5 than in the BQS, and fully compati- 
ble with other recent quasar surveys. 

Given apparent magnitude B j and redshift z, absolute mag- 
nitudes Mbj of the QSOs were estimated using the usual for- 
mula Mbj = Bj + 5 - 5\ogdL{z) +A+K{z), where dL{z) is 
the 'luminosity distance' in an expanding F riedm ann universe, 
computed with the relation given by Terrell ( 1977 ). We adopted 
Ho = 50kms^^ Mpc^^ and = 0.5 to specify the cosmologi- 
cal model. For the extinction term A, only attenuation by Galac- 
tic dust has been assumed, and the values of Abj were estimated 
from H I column densities as described in Paper 1 . 

Somewhat more complex is the situation for the K term: 
Although quasar spectra have often been approximated by sim- 
ple power laws /v °^ v" with adopted mean spectral index a 
between —1 and —0.3, actual spectral energy distributions are 
more complicated than this. We have recently determined a new 
K correction as a function of redshift that results in luminosi- 
ties lower by 0.4 mag for high-z QSOs, compared to previously 
published relations (Wisotzki 1999| ), and we adopt that new 
K{z) relation for the present paper The reader should be alert 
that the space densities derived in the next section are there- 
fore shifted towards lower luminosities, relative to what can be 
found in the literature. However, whenever comparisons with 
the results of earlier workers are carried out in this paper, in 
particular with parametric descriptions of the QLF and its evo- 
lution, we use the appropriate original K{z) relations employed 
by others. The qualitative conclusions of this paper do not de- 
pend on the details of the adopted K correction. 

Fig. |l] shows the distribution of sources over absolute mag- 
nitudes and redshifts. At z < 0.1, the majority of objects would 
usually be called Seyfert 1 galaxies rather than QSOs, but at all 
other redshifts there is no ambiguity of this sort, the HES al- 
ways sampling the optically brightest parts of the known quasar 
population. 



3. Luminosity functions 

Space densities at given absolute magnitude were computed 
in bins of redshift, and luminosity for the binned differential 
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Fig. 2. The cumulative quasar luminos- 
ity function determined from the HES 
sample for five adjacent redshift shells. 
Each panel shows the result for one red- 
shift domain, including error bars esti- 
mated from Poisson statistics. 



form, using the 1 /Vmax estimator (Felten 197^ . Five redshift 
intervals were chosen to represent the evolution of the bright 
end of the QLF with cosmic time, with boundaries z — 0.0-0.3, 
0.3-0.7, 0.7-1.25, 1.25-2.0, and 2.0-2.9. These intervals corre- 
spond to approximately equal steps in log(l +z) which is the 
cosmic time variable in the PLE model of Boyle et al. (1991), 



from where the above set of intervals was adopted. The lowest 
(z = 0-0.3) interval is not common in quasar evolution studies 
in the optical (see below), but has been added for this analysis. 

Fig. ^ displays the luminosity function, in cumulative form, 
for each shell separately. As discussed by Hewett et al. (1993), 
this way of representing a luminosity function has consider- 
able advantages over the conventional approach to bin the data 
in absolute magnitude: The contribution of each object is easily 
made apparent, the complicated weighting within the luminos- 
ity bins needs not be corrected for, and the full range of the data 
can be displayed without 'incomplete bin' edge effects. How- 
ever, the data are still binned in redshift and the constructed 
luminosity functions therefore do not represent the QLF at spe- 
cific epochs, but at 'mean' epochs that vary with luminosity 
(a form of Malmquist bias). In particular, in the presence of 
strong positive evolution, the apparent QLF for any given red- 
shift shell will appear considerably flatter than the intrinsic 



QLF. The statistical test procedures employed below are not 
affected by this bias. In the following we comment briefly on 
the observed properties of each of the subsets, 
z < 0.3: The local (z « 0) universe is the only domain where 
the full range of luminosities is technically accessible to a sin- 
gle flux-limited survey. The HES is the first optical quasar sur- 
vey capable of fully exploiting this option, being specifically 
designed to reduce selection biases due to the presence of ex- 
tended galaxy envelopes. This has already yielded a first es- 
timate of the local QLF based on a small subsample of the 
HES (Kohler et al. 1997). The present sample of 160 objects 



with z < 0.3 increases the number of 20 'local QSOs' used by 
Kohler et al. by almost an order of magnitude, allowing a much 
more accurate assessment of QLF shape and slope. The local 
luminosity function is discussed further in Sect, ^below. 
0.3 < z < 0.7: This shell is equal to the 'low-redshift' regime 
of most other optical quasar surveys. Host galaxies are visible 
only in exceptional cases, due to the (1 -l-z)'* surface brightness 
dimming and the steep K correction for galaxy spectra in the 
Bj band. For z < 0.5, quasar spectra are generally extremely 
blue and easily separated from their UV colours; around z — 
0.6, the short-wavelength end of the 'little blue bump' causes 
quasar colours to be somewhat redder, but it is shown in Paper 1 
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that there is no evidence for a significantly reduced degree of 
completeness in the HES sample in this redshift range. 
0.7 < z < 1.25: The redshift distribution seems to have a local 
maximum around z — 1.1, but this is unlikely to be a positive 
selection bias: There are no strong emission lines visible (only 
C III] A-1909 which is often just barely detected), and neither the 
colour dependence of z nor the K correction have discernible 
features at this redshift. The most probable interpretation is that 
of a statistical fluctuation. 

1.25 < z < 2.0: The bright end of the QLF in this shell is dom- 
inated by the exceptionally luminous QSO HE 0515—4415. In 
this redshift range, the total number of QSOs predicted by the 
B91 model is consistent with what is observed in the HES; 
however, there is a tendency that the number of fainter objects 
is predicted too high. This trend is continued, even stronger, in 
the LBQS as noted by Hewett et al. ( 1993 ), and it is also visible 
in the HBQS (La Franca & Cristiani 1997), and it is therefore 



likely that not incompleteness of the data, but an overprediction 
of the PLE model causes the discrepancy. 
2.0 < z < 2.9: The high-redshift range of the present HES 
sample is still only sparsely populated. Up to z ~ 2.6, the 
HES selection criteria are predominantly sensitive to contin- 
uum properties. Although Lya is generally detected as a strong 
emission line from z ^ 1.8 onward, the broad range of colour 
criteria employed in the HES ensures that all of the z < 2.6 
QSOs in the sample are in fact colour-selected (cf. Paper 1), 
thus avoiding the strong biases associated with emission line 
detection. The expected quasar colours as functions of z show 
a strong maximum of UV excess around z — 2.\- therefore the 
decrease of numbers when passing the z = 2 mark cannot be 
related to selection effects. 

At z > 2.6, the QSOs were generally found by feature- 
detection algorithms, depending on the presence of a strong 
La emission and/or a pronounced spectral break at the onset 
of the Lyman forest. Both features are common but not neces- 
sarily strong, making it hard to quantitatively assess the degree 
of (in)completeness. Another difficulty is the increasingly ill- 
defined K correction, with intergalactic Lya forest absorption 
adding to the uncertainty of spectral energy distribution. The 
number of objects in the sample should be, at any rate, con- 
sidered as a lower limit to the number expected for a 'perfect' 
sample limited only by By magnitude. 

In Fig. ^ all data have been combined into one frame, pre- 
sented in the traditional binned differential form. The HES pro- 
vides substantial improvement over previous surveys in cov- 
ering the bright end of the QLF at all redshifts z ^ 3, plus a 
full construction of the important local luminosity function. 
For comparison, the dotted line in Fig. |] shows the differen- 
tial QLF, integrated over the appropriate redshift ranges, from 
the Boyle et al. (|1991|) PLE model. 



4. Constraints on parametric models 

4.1. Pure luminosity evolution 

By combining the ^ 400 QSOs from the AAT multifibre sur- 
vey with a set of other, smaller quasar samples - in particu- 




Fig. 3. Evolution of the differential QLF in bins of 1 mag with 
redshift. The redshift shells are the same as in Fig. ^ from left 
to right. The dotted line shows the prediction of the Boyle et al. 



(1991) PLE model. 



lar the BQS at the bright end -, Boyle et al. ( |1988| ) suggested 
that the quasar luminosity function can be approximated as 
a (logarithmically) shape-invariant double power law, merely 
shifted in luminosity by an offset governed through a redshift- 
independent evolution parameter ki w 3-3.5. In a later analysis, 
Boyle et al. (1991; hereafter B91) substantiated this claim and 
extended the PLE model up to z = 2.9, albeit with the introduc- 
tion of an additional parameter Zmax above which the comoving 
space densities were kept constant. 

Comparing the prediction of the specific parametric model 
of B91 with the empirical LF estimates based on the new HES 
sample (Fig. ||) reveals substantial discrepancies at most red- 
shifts: 

- At low and intermediate z, there is a significant excess of 
high-luminosity QSOs over the PLE prediction. 

- At low redshifts, there are less intermediate-luminosity 
QSOs than demanded by PLE. 

- The characteristic 'break' observed at high z does not ap- 
pear in the local (z < 0.3) QLF. 

- The QLF is much steeper at high redshift. 

- There is evidence for continued evolution after z > 2. 

Most of these discrepancies, in particular the feature of a 
relatively flat low-redshift QLF, were already detected in pre- 



vious surveys (Hewett et al. 1993; Miller et al. 1993) and are 



clearly not artefacts of the HES selection procedure; the HES 
confirms these findings with high significance. Because the 
actual low-z QLF is flatter than the PLE prediction, discrep- 
ancies show up most prominently at the highest luminosities. 
With its large effective area and the luminosity range probed, 
the HES provides considerably improved statistical coverage to 
test the validity of the PLE parametrisation. A two-dimensional 
Kolmogorov-Smirnov test (Peacock 1983 ) applied to the joint 
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In a recent attempt that includes several sets of new survey 
material come available in the 1990s, La Franca & Cristiani 
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Fig. 4. Variation of the bright-end slope yi with redshift. 



distribution of (z, M) pairs yields a formal acceptance proba- 
bihty oiP = 10"^ for the redshift range 0.3 < z < 2.2, and P = 
10"'^ for the whole definition interval of the HES, < z < 3.2, 
so the B91 model is rejected with very high confidence. 

An essential feature of all PLE models is the fact that the 
slope at the bright end of the QLF, yi , is necessarily constant. 
Fig. ^ shows that except for the z < 0.3 shell, the HES objects 
are always located on the steep part of the QLF, well beyond 
the ubiquitous 'break' ; this will be so for all realistic PLE mod- 
els. A good test for the validity of PLE models in general can 
be made by monitoring the bright-end slope of the QLF as a 
function of z. To this purpose we have fitted straight lines to 
the cumulative logarithmic LF within each redshift shell (i.e., 
assuming a power law). The result is shown in Fig. ^ where 
yi(z) is resolved into redshift bins of Alog(l +z) = 0.05. Al- 
though Malmquist bias and evolution within the shells tend to 
make the measured slope at given z somewhat flatter than the 
intrinsic one, the bins are small enough to ensure that this bias 
is almost negligible. 

The trend for the bright-end slope to become flatter with 
decreasing redshift is very obvious, confirming the earlier con- 
jecture by Goldschmidt & Miller ( |1998 ), with significantly im- 
proved redshift resolution. Fig. ^ also shows that this trend is 
confined to the redshift regime z ^ 1 , while for z > 1 the slope 
is consistent with being constant at yi ~ —3.7, approximately 
equal to the va lue of yi in most PLE fits (B91; La Franca & 
Cristiani 1997 1. Fig. Qmay be taken to furthermore suggest that 
the QLF steepens again after z > 2, but that is not statistically 
significant. 

This test demonstrates that the basic assumption of all PLE 
models - a redshift-invariant QLF shape - is in conflict with ob- 
servations, independently of the details of the parametrisation. 
The bright end of the quasar luminosity function undergoes a 
gradual transformation, becoming flatter with cosmic time. 



(1997) modelled the QLF with a 'luminosity-dependent lumi- 
nosity evolution' (LDLE) scheme, by equipping the PLE for- 
mula of Boyle et al. with an additional term to produce a flat- 
ter QLF at low redshift while maintaining the PLE description 
for higher redshifts. We have subjected their model (includ- 
ing the correction from La Franca & Cristiani 1998) to similar 
tests as the Boyle et al. model, with the following results: (1) 
There is very little difference in performance between LDLE 
and PLE for z ^ 1 . (2) The bright end slope yi decreases (by 
absolute value) towards lower redshifts as demanded by the ob- 
servations, but the number of luminous QSOs is overpredicted. 
(3) Below z — 0.4 the LDLE formalism assigns strongly nega- 
tive evolution to the most luminous QSOs, producing an almost 
flat (yi = — 1) local luminosity function, which is of course in- 
compatible with all available data. In conclusion, the LDLE 
scheme as proposed by La Franca & Cristiani does not resolve 
the shortcomings of PLE, and it is completely inadequate for 
the low-redshift domain. 



4.3. Pure density evolution 

Although pure density evolution (PDF) is the oldest and intu- 
itively simplest form to parameterise an evolving luminosity 
function, it has been ruled out by several independent investi- 
gations. The argument is usually based on the counts of faint 
QSOs at high redshifts, which are observed to be much rarer 
than predicted by PDF. Even if faint QSO surveys are assumed 
to be substantially incomplete, a stringent limit is set by the 
measured intensity of th e diffu se extragalactic soft X-ray back- 
ground (Marshall et al. 1983 ). In a preliminary ana lysis of a 
composite ROSAT-selected AGN sample, Hasinger ( 1998 ) re- 
cently challenged the unequivocal rejection of the PDF con- 
cept, claiming that simple density evolution in fact provided a 
much better fit to the data than PLE. This claim has meanwhile 



been qualified (Miyaji et al. 1998 ; see below), but the perfor- 
mance of PDF with these new samples seems indeed surpris- 
ingly good. 

While a detailed discussion of the similarities and differ- 
ences between optically and X-ray selected samples is clearly 
outside the scope of this paper, it may be worthwhile to in- 
vestigate the merits of PDF for the bright HES sample. Fit- 
ting a one-parameter PDF model to the data is simple: We use 
the modified variable V'/^max defined as the density-weighted 
volume integrals, V' = / p(z)t/V, where p(z) is the density 
evolution function. Adopting the customary form for p(z) = 
p(0) (1 +zY", ko is varied until the sample average of y'/y^^^ 
reaches 1/2. Based on the new K correction and using the full 
HES sample within < z < 3.2, this yields ko = 5.35, thus a 
much smaller value than obtained in earlier fits based on fainter 



samples, but very similar to the X-ray value of Hasinger (1998 ). 
The distribution of y '/^max is also consistent with being uni- 
form within [0,1]. However, this is only a necessary and not a 
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sufficient condition; the two-dimensional KS statistic yields a 
formal acceptance probability of P = 10^'^. For the first time, 
PDE can be significantly rejected from a bright QSO sample 
alone, without any reference to faint QSO counts. The reason 
is, of course, again the steepening of the bright QLF part to- 
wards higher redshift, leading to incompatible QLF shapes at 
low and high z- Note that although the z — Q and the z 2 QLF 
sections as seen in the HES do not overlap in luminosity, there 
is sufficient dynamic range at intermediate redshifts to monitor 
the steepening slope at given Mbj ■ 

4.4. Luminosity-dependent density evolution 

A luminosity function changing with cosmic time can be 
always be described by the general concept of luminosity- 
dependent density evolution (LDDE). Concentrating on the 
more luminous part of the QSO population, Schmidt & Green 
( 1983 ) proposed a specific LDDE scheme in which the more lu- 
minous QSOs undergo stronger evolution. (Note that this is not 
necessarily in conflict with the principles of PLE, which could 
be translated into an equivalent LDDE formalism.) However, 
in the Schmidt & Green picture, the QLF becomes successively 
steeper at smaller redshift, opposite to what is observed in the 
new samples. 



A related approach was followed by Miyaji et al. ( 1998 ), 
accommodating the proposed density evolution of soft X-ray 
selected AGN (Hasinger 1998 ) with the constraints set by the 
X-ray background. In their formulation, the density evolution 
index ko{L) is smallest on the low-luminosity tail of the QLF; 
above a certain luminosity threshold, the index remains con- 
stant. The latter feature implies that directly transfering this 
model to the evolution of optically selected QSOs would re- 
sult in the same conflicts as for the case of PDE discussed in 
the preceeding subsection. The specific merits of the Miyaji et 
al. LDDE model are in the low-luminosity regime, where the 
HES yields no constraints on the evolution rates, apart from the 
local luminosity function. 

While a single survey such as the HES provides insuffi- 
cient coverage of the Hubble diagram to allow new parametric 
models to be developed, it seems safe to conclude from the 
present data that only a rather complex luminosity-dependent 
evolution scheme will be capable to reproduce all observed 
features, making quasar evolution more complicated than pre- 
viously thought. 

5. The local luminosity function 

The basic results of our earlier analysis of the local QLF 
(Kohler et al. 



1997 



hereafter K97) were: (1) The number of 
high-luminosity low-redshift QSOs has previously been under- 
estimated. (2) The local QLF does not display any evidence of 
a break or significant change of slope at absolute magnitudes 
Mbj ^ —20. The new HES sample improves the K97 data by a 
factor of 6 in covered area and by a factor of 8 in sample size. 
Both above mentioned results are confirmed in the present pa- 
per, although the discrepancy in number counts is less extreme 



than found in the K97 sample. However, considering the error 
bars, the results of the two analyses are fully compatible. 

A new feature needing corroboration is the apparent flatten- 
ing of the local QLF at luminosities fainter than Mgj — —20, 
apparent in Fig. ^ Note that it is based on only 3 sources 
with Mbj > —20, thus the statistical significance of the fea- 
ture is very low. We are currently investigating this part of the 
local QLF in more detail, including a dedicated treatment of 
host galaxy influences and including the transition to low-level 
AGN in nearby galaxies. For the present discussion we concen- 
trate on objects with Mbj < —20. Even for these, host galaxy 
contributions to standard isophotal magnitudes would be non- 
negligible. To avoid this effect, the HES magnitudes were mea- 
sured over effective apertures of the size of the seeing disk, so 
the absolute magnitudes do not correspond to total, but effec- 
tively to nuclear luminosities; even large and bright host galax- 
ies do not contribute with more than their seeing-convolved 
central surface brightness to the measured flux. Thus, while the 
measured magnitudes are not individually corrected for host 
galaxy contributions as it was the case in K97, they are never- 
theless dominated by the active nuclei, not by the galaxy hosts. 

The accurate determination of the local luminosity function 
of QSOs and Seyfert nuclei is an important, but also a difficult 
task. In addition to the challenge of generating unbiased sam- 
ples, there are several problems to be overcome: Disentangling 
the contributions of nuclei and hosts to the measured luminos- 
ity is a problem not only in the optical domain, but also at X- 
ray wavelengths (Lehmann et al. 1998 ). A further complication 
arises from the probable bifurcation of the AGN population 
into a Seyfert 1 and a Seyfert 2 branch. These points need to 
be addressed before a meaningful physical interpretation of the 
local QLF can be made. 



6. The most luminous QSOs at high redshift 

6.1. Evolution of space densities 

The substantial uncertaintes in the present knowledge of quasar 
evolution at redshifts z > 2 are mostly related to the increas- 
ing difficulties of obtaining sizeable well-defined flux-limited 
samples. The traditional UV excess method based on U — B 
colours breaks down at z = 2.2 where the Lya emission line 
is redshifted into the B band; multicolour techniques are ca- 
pable of reaching higher z, but with the penalty of a compli- 
cated redshift-dependent selection function. Photographic ob- 
jective prism surveys such as the LBQS and the HES, on the 
other hand, perform quite well up to z ~ 3, but become rapidly 
incomplete towards z = 3.4 where Lya moves out of the de- 
tector bandpass (see also Sect. ||). While the usefulness of the 
HES sample for space density estimation is therefore limited 
to z ^ 3, this is just the redshift region where most previous 
workers located the maximum of quasar activity. Furthermore, 
multicolour surveys show greatly reduced performance in this 
z range, including the forthcoming SDSS (Fan 1999). 



While at low redshift the sampled luminosity range of a 
flux-limited survey depends sensitively on z, this dependency 
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Fig. 5. Evolution of the integrated space density for QSOs with 
Mb < —28. The filled symbols show the values computed from 
the HES, the dotted line delineates the prediction of the B91 
PLE model. Open symbold indicate the available higher red- 
shift measurements for this luminosit y reg ime (triangle: from 
Hewett et al. 



et al. 



1995) 



1993; circle: Irwin et al. 1991; square: Kennefick 



is much weaker towards higher redshift, and here it is possible 
to trace quasar space densities/or given luminosity over a wide 
redshift range. The HES provides appropriate data material to 
apply this to the rare species of very high luminosity QSOs, 
following their evolution up to z ~ 3. 

Fig. H shows the evolution of integrated comoving space 
density for Mbj < —28 over several disjoint redshift bins. Be- 
yond z > 2.5, even the faintest HES quasars are more luminous 
than this. The last plotted value is therefore based on a linear 
fit to the (logarithmic) cumulative luminosity function in the 
corresponding bin, extrapolated to Mbj = —28. This increases 
the uncertainty, but does not bias the estimate. Recall also that 
QSO selection gets increasingly difficult close to the HES high- 
redshift limit, so that the true space density for this last point 
will be even higher. 

The steady increase in space densities, up to the highest red- 
shifts sampled, is remarkable. In particular, there is no indica- 
tion for a turnover or only a slowing down of evolution. Unless 
the space density as a function of redshift shows an extremely 
narrow peak just below z = 3 - which would still be consis- 
tent with our data, considering the error bars -, we conclude 
that the maximum space density has almost certainly not yet 
been reached by our survey and must be located well beyond 
z = 3. This is in marked contrast to the general notion that the 
maximum of quasar activity occurs at z < 3 (e.g., Schmidt et 
al. 1995; Shaver et al. 1996). Two possible interpretations are 



conceivable: 

(1) The maximum of quasar activity might depend on lu- 
minosity, in the sense that high-luminosity objects have a max- 
imum shifted towards higher redshifts. Only three published 



surveys can provide meaningful additional datapoints for Fig. 
5|: Th e estimate for 3 < z < 3.4 from the LBQS (Hewett et al. 
1993) which is almost certainly only a lower limit; the z > 4 



survey from Irwin et al. (1991), and the similar sample from 



Kennefick et al. (1995). None is in direct conflict with the HES 



data, but the high space density claimed by Irwin et al. (1991) 



seems to be joining particularly well with the HES. However, 
the origin of the discrepancy between the results of Irwin et al. 
and Kennefick et al. - based on partly the same objects - is not 
sufficiently understood. Kennefick et al. suggested that differ- 
ent adopted K corrections may be responsible, but this has not 
been confirmed so far 

Additional support for a stronger evolution of high- 
luminosity QSOs comes from the analysis of Warren et al. 



(1994): although their survey does not reach the luminosity 



range in question, their adopted luminosity evolution model in 
which only the bright part evolves can be extrapolated to evo- 
lution rates around z — 2-3 that are quite compatible with the 
HES data. 

If it should prove correct that high-luminosity QSOs con- 
tinue to show positive evolution where the space densities of 
their lower luminosity counterparts already turn over, this has 
substantial consequences for cosmogony and galaxy formation 
scenarios. High luminosity presumably implies host galaxy 
mass (in analogy to the confirmed such relation at low redshifts. 



cf. McLeod & Rieke 1995). This would impose significant con- 



straints on theories for the formation of very massive structures 
in the early universe. 

(2) Incompleteness of fainter QSO surveys provides an 
obvious alternative explanation. Several recent deep surveys 
found high-redshift QSOs in larger numbers than anticipated 
(Miyaji et al. |1998|; Wolf et al. |1999[), although intercompar- 



isons are difficult as each survey samples a different luminosity 
regime, and the statistics are still very poor. 

6.2. Gravitational lensing 

Statistical analyses based on optical quasar surveys can suf- 
fer from biases quite similar to those present in low-redshift 
QSO host galaxy studies: Systematic rejection of non-pointlike 
sources, a property of most surveys, invariably leads to miss- 
ing large-separation lenses. The HES avoids this selection bias, 
and since since it samples the region of the QLF where 'mag- 



nification bias' (Turner 1980) is expected to be most effective. 



it is well suited as a testing ground for lensing statistics. While 
a full investigation of this subject is beyond the scope of this 
paper, we can make a few qualitative inferences on the effects 
of lensing onto the observed luminosity function. 

The principal result of 'magnification bias' is to create a 
flatter observed QLF, compared to the intrinsic one. It can be 
shown that the effect is noticeable only when the QLF slope 



is steep (e.g., Schneider et al. 1992 1. As the optical depth for 
gravitational lensing increases with z, the magnification bias 
is expected to be stronger for sources at higher redshift. Un- 
der the assumptions of an intrinsically constant QLF slope and 
magnification bias being a dominant effect, we should observe 
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a bright end of the QLF that flattens with increasing redshift, 
opposite to what is actually observed. This does not exclude 
that there might be some influence of lensing on the observed 
luminosity function, but the observations constrain magnifica- 
tion bias to be of very minor importance for the QLF as a whole 
up to z ~ 2. This is consistent with the fact that only ~ 1 out of 
100 QSOs with My < —28 shows discernible image splitting 



due to lensing (Surdej et al. (1993) 



7. Conclusions 

The analysis of a new sample of very bright QSOs has shown 
that quasar evolution is more complex than generally as- 
sumed. The Hamburg/ESO survey provides insight into oth- 
erwise rarely sampled regions of the luminosity-redshift plane, 
enabling a direct assessment of the local luminosity function, 
and to probe the evolution of the bright part of the QLF up 
to z ~ 3. While the high-redshift QLF has a very steep high- 
luminosity tail, the low-redshift QLF approaches towards a 
simple power-law form with rather flat slope. No simple 'pure' 
evolution scheme based on a shape-invariant QLF is capable of 
reproducing this transformation. Future analyses of large com- 
posite samples covering a wide range of luminosities are re- 
quired to show if appropriate parametric descriptions can be 
found that are consistent with this behaviour. 
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